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Abstract. We introduce a general-purpose framework for interconnecting scientific 
simulation programs using a homogeneous, unified software interface. Our framework 
is intrinsically parallel, and conveniently separates all components in memory. It per- 
forms unit conversion between different modules automatically and defines common 
data structures to communicate across different codes. We use the framework to sim- 
ulate embedded star clusters. For this purpose we couple solvers for gravitational dy- 
namics, stellar evolution and hydrodynamics to self consistently resolve the dynamical 
evolution simultaneousy with the internal nuclear evolution of the stars and the hydro- 
dynamic response of the gas. We find, in contrast to earlier studies, that the survival 
of a young star cluster depends only weakly on the efficiency of star formation. The 
main reason for this weak dependency is the asymmetric expulsion of the embedding 
gas from the cluster. 



1. Introduction 

Large-scale, high-resolution computer simulations dominate many areas of theoreti- 
cal and computational astrophysics. The demand for such simulations has expanded 
steadily over the past decade, and is likely to continue to grow in coming years due 
to the relentless increase in the volume, precision, and dynamic range of experimen- 
tal data, as well as the ever-widening spectral coverage of observations. In order to 
accommodate the improved observations, simulations must see a comparable improve- 
ment in detail. In recent years, simulation environments have grown substantially by 
incorporating more detailed descriptions of more physical processes, but the fundamen- 
tal design of the underlying codes has remained un changed sin ce the introduction of 
object -oriented programming iMcCarthy et all (fl 962) and patterns [Kent & Cun ningham 



(119871) . As a result, maintaining and extending existing large-scale, multi-physics 
solvers has become a major undertaking. The legacy of design choices made long ago 
can hamper further code development and expansion, prevent scaling on large parallel 
computers, and render maintenance almost impossible. 

The root cause of the increase in code complexity lies in the traditional approach of 
incorporating multi-physics components into a single simulation — namely, solving the 
equations appropriate to all components in a monolithic software suite, often written 
by a single researcher. This monolithic solution may seem desirable from the stand- 
point of consistency and performance, but the resulting software generally suffers from 
fundamental problems in maintenance and expansion. 
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The AMUSE (Astrophysical Multipurpose Software Environment^] project as- 
similates well tested applications into a software suite with which we can perform in- 
dividual tasks, or reassemble the parts into a new application that combines a wide 
variety of solvers. The interfaces of codes within a common domain are designed to be 
as homogeneous as possible. This approach is possible in astrophysics because of the 
tradition among astronomers during the last several decades of sharing scientific soft- 
ware. Many of these applications were written by experts who spent careers developing 
these codes and using them to conduct a wide range of numerical experiments. These 
packages are generally developed and maintained independently of one another. We 
refer to them collectively as "community" software. AMUSE has rece ntly surpassed 
our "Noah's ark" developmental milestone (Porteg ies Zwart et al. 2009), in which we 
have at least two numerical solvers for each of the astrophysical domains of interest: 
gravitational dynamics, stellar evolution, hydrodynamics, and radiative transfer. 

In a first step towards using this framework, we combine here three of these fun- 
damental ingredients of AMUSE to address a long-standing problem in astrophysics: 
the relevance of the star formation efficiency to the survival of embedded star clusters. 



2. AMUSE 

The AMUSE environment allows astrophysical codes from different domains to be 
combined to conduct numerical experiments. The community codes are generally writ- 
ten independently, so AMUSE encompasses a wide variety of computer languages and 
programming styles. The fundamental design feature of the framework is the abstrac- 
tion of the functionality of individual community codes behind physically motivated in- 
terfaces that hide their complexity and/or numerical implementation. AMUSE presents 
the user with standard building blocks that can be combined into applications and nu- 
merical experiments. 

The binding language that stitches the codes together is Python. The relatively 
low speed of this high-level language is not an issue, since the focus in the high-level 
management code is not so much performance (the computational cost being concen- 
trated in the component codes), but algorithmic flexibility and ease of programming to 
allow rapid prototyping. As described in more detail in the contribution by McMillan, 
Portegies Zwart, and van Elteren elsewhere in these proceedings, an AMUSE applica- 
tion consists of a Python user script controlling one or more community modules. The 
user script specifies the initial conditions, selects the simulation modules, and manages 
their use. The coupling between the user script and a community code is handled by 
the community module, which contains an MPI-based communication interface onto 
the code, as well as unit-handling facilities and an object oriented data model. 

The relationships among the community codes define the numerical experiment. 
Our model here combines the effects of the self-gravity and nuclear evolution of the 
stars with the hydrodynamics of the intracluster gas (Pelupessy & Portegies Zwart 201 1, 
in preparation). The latter includes both the primordial gas content of the cluster and 
the gas liberated by the stars via stellar winds and supernovae. In this case we construct 
a hybrid N-body/stellar/hydrodynamic solver by combining a direct N-body integrator, 
a stellar evolution package, and an SPH code. 



see www . amusecdoe . org. 
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The Python user script controlling the experiment generates the initial conditions 
(masses, positions and velocities of the stars, and the distribution of the gas), specifies 
the various solvers, structures the procedural calling sequences, resolves all interactions 
among the various physical domains (e.g. feedback from the stellar winds and super- 
novae to the surrounding gas), and processes the out put. The particu lar modules em- 
ployed in this experi ment are the Gadget -2 SPH code ISpring"ell (120051) . t he PhiGRAPE 
Hermite N-body code lHarfst et al.1 (120071), the tree gravity code Octarav lGaburov et al. 



2010) and the stellar evolution code iHurley et al.l d2000). The first dynamical model 
is used for the integration of the equations of motion of the stars, the second mod- 
ule is used for the gravitational coupling between the gas-particles and the stars. The 
combined solver consists of an integrator for the coupled gas/gravitational dynamics 
systems and a feedback prescription for mechanical energy input from the evolving 
stars. 

The gas and gravitational dynamics are coupled via the BRIDGE integrator (iFujii et al. 



2007). BRIDGE provides a semi-symplectic mapping for gravitational evolution in 



cases where the dynamics of a system can be split into two (or more) distinct regimes. 
A typical application would be a dense star cluster in a galaxy, where the internal dy- 
namics of the former evolves on a relatively shor t timescale compared to t he dynamics 
of the latter. A similar idea was implemented by ISaitoh & Makinol (1201 Oh by splitting 
the gravitational and hydrodynamic evolution operators for simulating gas-rich galaxy 
mergers. They expressed the algorithm in a single monolithic code, whereas we adopt 
the concept of operator splitting within AMUSE to couple different codes. 



3. Initial conditions 



The clusters we simulate are co mposed of a mixture of gas and N = 1000 stars; both are 
distributed in a lPlummerl (Il91lh s phere, and they have the same characteristic radius. 
Stellar masses are assigned using a lSalpetej (Il955h IMF between 0.1 and 100 M , with 



an additional constraint that the most massive star is ~ 22 M©. This maximum mass is 
based on the most massive star naively expect ed for a cluster with this number of stars 



and mass function iKroupa & Weidnen (l2003h . The masses of the stars are assigned 



independently of their positions in the cluster. We present here the results of two of our 
simulations, which in our larger paper describing this work are identified as model A2 
and model A5 (see Pelupessy & Portegies Zwart 2011). 

For small clusters the number of high-mass stars can vary quite substantially be- 
tween different realizations of the IMF, and we have performed simulations of models 
A2 and A5 with numerous random realizations of the IMF to examine this effect. We 
have performed additional simulations in which we varied the number of gas particles 
to test whether our results are independent of the resolution of the gas dynamics. 



4. Results 

FigureQ] shows the stellar and gas distribution of our models A2 and run A5. In both 
models we parameterize the relative feedback efficiency between the stars and the gas 
by a parameter /ft>, which is the fraction of the total supernova and wind energy output 
that ends up as thermal energy in the ISM (this accounts for the uncertainties in mod- 
eling the feedback and radiative losses). For the A2 model we take = 0.1 while 
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for the A5 model fa = 0.01. The feedback is implemented by returning gas parti- 
cles from stars in proportion to the mass loss rates of the stellar wind and SN, with a 
thermal energy set by the mechanical luminosity of the star and the canonical energy 
E sn = 10 51 ergs in case a star goes supernova. For each model we plot the stars and a 
slice through the gas density distribution at four moments during the simulation. In the 
first A2 frame (at 0.96 Myr), we see the early stages where stellar winds create buoyant 
bubbles that rise out of the potential of the star cluster. As the mechanical luminosity 
increases these bubbles grow until they blow away sizable fractions of the cluster gas 
and a free-flowing wind develops (4.37 Myr frame). The strong feedback then unbinds 
most of the gas of the cluster. At approximately 9.5 Myr the cluster ISM has been 
ejected — the gas visible in this frame originates from the strong AGB wind of the most 
massive progenitor (m ~ 21 M Q ). 

At an age of 9.54 Myr the most massive star in the simulation undergoes a super- 
nova explosion which ejects the remaining gas from the cluster (both simulations use 
the same initial realizations for the IMF and stellar positions). For the A5 simulation 
(with a relative feedback efficiency of 0.01, compared to 0.1 for model A2), the initial 
wind stages (before 0.96 Myr) proceed less violently, with smaller bubbles, and a free- 
flowing wind does not develop until just before the supernova (compare the 4.37 Myr 
frames). The main difference between the A2 and A5 runs is that most of the cluster 
gas is retained in the latter case until the first supernova blows it away. Just before the 
supernova the A5 cluster is much more compact than in the A2 run. 

Loosely bound associations can be distinguished from 'true' stellar clusters using 
limited observables by conside ring the ratio II of age T c \ and the crossing time T cross 



Gieles & Portegies Zwartl (1201 lh : 

n = r c i/r cross (i) 

In Fig. [2] we present the value of II for models A2 and A5 as functions of time. A value 
of II < 1 indicates that the cluster is in a ballistic state of expansion, whereas a value of 
II > 1 implies a bound state. In Fig.[2]we see two completely different behaviors for the 
evolution of IT. For model A2, II rises sharply but the cluser fails to reach a bound state 
(in the sense of II > 1) until much later (after t > 30 Myr), whereas model A5 reached 
a bound state within a few Myr after formation and remained marginally bound for the 
stof the simulation (up to ~ 30 Myr). In our survey of parameter space (see Pelupessy 
& Portegies Zwart 2011), we explore a wider range of initial conditions. 



5. Conclusions 

We have simulated star clusters in their embedded phase. Our simulations include the 
gravitational dynamics of the stars, the dynamics of the intracluster gas, and the internal 
evolution of the stars. We find that the star formation efficiency is a poor predictor of 
final state of the cluster. There are several arguments why the star formation efficiency 
is less important than has been found in earlier studies. The most dramatic event in the 
lifetime of a young cluster is the occurrence of the first supernova, which blows away 
most of the residual gas in the cluster. But due to earlier fast Wolf-Rayet winds from the 
massive stars most of the gas has already escaped without much damage to the cluster. 
In addition, during the time between the strong Wolf-Rayet wind and the supernova 
explosion the cluster has time to relax, making it more resilient against destruction by 
the loss of primordial gas. 




Figure I. Stellar and gas distribution of the A2 and A5 runs. The left panels 
show the gas and stellar distribution of the A2 run, those on the right panels the A5 
run. Snapshots are labeled by time in the lower right corner. The density plots show 
cuts through the mid-plane. The points show stars in 4 mass groups (m < 0.9 M©, 
0.9M o < m < 2.5M Q , 2.5M <m< 10M o andm> 10M o ). 
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time time 



Figu re 2. The ratio of the cluste r age to the calculated crossing time as defined 
by iGieles & Portegies Zwartl (1201 lb (II = age/fhc) f° r models A2 and A5 (see also 
Fig.©. 



Statistical variations in our method of generating the initial mass function have a 
profound effect on the early evolution of the cluster. The survival of the star cluster 
may well depend on the masses and orbits of the few most massive stars it contains. A 
cluster with a slight enhancement of massive stars may well dissolve, whereas a more 
fortunate cluster may be born with a larger gap between the masses of few most massive 
stars. Slight differences of even a few M Q in the most massive stars may well be crucial 
in determining the survival of the cluster. 

The surviving clusters are strongly mass segregated. During the embedded phase 
massive stars easily sink to the cluster center. The degree of mass segregation found in 
the surviving clusters nicely matches those required to explain the oberved degree of 
mass segregation in the Pleiades. 

Our prescription for the radiative feedback in our models is still very limited, and 
the next obvious step in improving our model would be by adopting a radiative transfer 
code to resolve this problem. 

Acknowledgments 

This work was supported by NWO (grants #643.200.503, #639.073.803 and #614.061.608), 
NOVA and the LKBF in the Netherlands, and by NSF grant AST-0708299 in the U.S. 



References 

Fujii, M., Iwasawa, M., Funato, Y, & Makino, J. 2007, PASJ, 59, 1095. 0706 . 2059 

Gaburov, E., Bedorf, J., & Portegies Zwart, S. 2010, ArXiv e-prints. 1005 .53841 

Gieles, M., & Portegies Zwart, S. F. 2011, MNRAS, 410, L6. I1010 1720 

Harfst, S., Gualandris, A., Merritt, D., Spurzem, R., Portegies Zwart, S., & Berczik, P. 2007, 

New Astronomy, 12, 357. |arXiv: astro-ph/0608125| 

Hurley, J. R., Pols, O. R., & Tout, C. A. 2000, MNRAS, 315, 543. |astro- ph/0001295 
Kent, B., & Cunningham, W. 1987, in OOPSLA '87 workshop on Specification and Design for 

Object-Oriented Programming. URL h ttp : //c2 . com/doc/ oopsla87 . html| 
Kroupa, P., & Weidner, C. 2003, ApJ, 598, 1076 



7 



McCarthy, J., Abrahams, P., Edwards, D., Hart, S., & Levin, M. I. 1962, LISP 1.5 Pro- 
grammer's Manual: Object - a synonym for atomic symbol (MIT Press), 105. URL 

http : //community . computerhi story . org/scc/pro j ects/LISP/book/LISP%2® 1 . 5%2 ©Programmer s%2®Manu 
Plummer, H. C. 1911, MNRAS, 71, 460 

Portegies Zwart, S., McMillan, S., Harfst, S., Groen, D., Fujii, M., Nuallain, B. 6., Glebbeek, 

E., Heggie, D., Lombardi, J., Hut, P., Angelou, V., Banerjee, S., Belkus, H., Fragos, T., 

Fregeau, J., Gaburov, E., Izzard, R., Juric, M., Justham, S., Sottoriva, A., Teuben, P., 

van Bever, J., Yaron, O., & Zemp, M. 2009, New Astronomy, 14, 369. 0807 . 1996 
Saitoh, T. R., & Makino, J. 2010, PAS J, 62, 301. 1Q9Q8 . 1 46Q 

Salpeter, E. E. 1955, ApJ, 121, 161 

Springel, V. 2005, MNRAS, 364, 1105. |arXiv: astro -ph/Q5Q5Q!Q 



user script 



1 rs 



fcuO 
O 




-0.5 



0.0 



0.5 1.0 

$log(t/Myr) 



1.5 



2.0 



community code 



Send ar swer 



3St 



2.5 



orbital_perio 



